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I. INTRODUCTION 



A. TARGET IDENTIFICATION 

In most radar applications so far, the only properties 
of the target that are measured are its location in range 
and angle. Such radars are sometimes called blob detectors 
since they recognize targets only as "blobs" located some- 
where in the space. In the era of modern electronic warfare 
it is desirable, however, to extract more information about 
the target. Radar should be able to recognize one type of 
target from another; that is, to determine that the target is 
a Boeing 747 aircraft and not a DC- 10, or that the partic- 
ular ship is a warship and not a freighter. This capability 
is called target identification or classification. 

Generally, radar target identification requires a larger 
signal to noise ratio then normally needed for detection. 
Hence, it is to be expected that the range at which target 
identification can be made would be less than the range at 
which the target can be first detected. 

One of the approaches that might be used for target 
identification is that of electromagnetic back scattering. 

The frequencies that are suggested as being important 
are those corresponding to wavelengths from half the target 
size to wavelengths about one-tenth the target dimensions 
The use of such frequencies is said to provide overall 
dimensions, approximate shape, and material composition of 
the target [Ref. 1]. 

It is particularly useful to investigate the resonances 
excited by electromagnetic fields of different frequencies 
incident on the target. At resonant frequencies the scat- 
tered fields are stronger compared to non-resonant 
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situation. The back scattering cross sections are larger 
and provide easier target detection and identification. 

Associated with resonant frequencies are the natural 
modes. A natural mode is defined as a possible "source- 
free" distribution of surface current on the scatterer 
according to Stratton [Ref. 6]. When the target is irradi- 
ated with an incident electromagnetic wave, surface current 
is excited on it. This surface current radiates and gener- 
ates the scattered field. A resonance occurs when one or 
certain combination of the natural modes responds very 
strongly to a particular type of excitation. The natural 
modes are defined only by the geometrical structure and the 
composition of the target, hence the resonant frequencies 
will be only functions of the target's physical configura- 
tion. Note that the resonant frequencies refered to in this 
thesis are real constants, in contrast to the "natural reso- 
nance" frequencies of the Singularity Expansion Method 
[Ref. 2]. 

Investigation of target back scattering over a broad 
frequency band appears to be a promising method to achieve 
target identification. If the phase shift as well as the 
amplitude of the scattered field are measured, fewer 
frequencies might be utilized than when amplitude alone is 
obtained [Ref. 3]. 

The polarization of the radar back scattered energy 
depends on the target properties and differs, in general, 
with the polarization of the energy incident on the target. 
This property can be used as a possible basis for discrimi- 
nating one target from another [Ref. 4]. For example, a 
thin straight wire can be readily distinguished from a homo- 
geneous sphere by observing the variation of the back scat- 
tered field amplitude as the polarization is rotated. The 
field amplitude from the sphere will be unmodulated, and the 
back scattered field amplitude from the wire will vary 
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between a maximum and minimum at twice the rate at which the 
polarization is rotated. By transmitting two orthogonal 
polarizations and measuring the amplitude and phase of the 
received back scattered field on each polarization, as well 
as the cross-polarization component, additional means of 
target identification can be provided. It is worthwhile to 
indicate that non-cooperative target recognition through 
various methods has been an ongoing topic of intensive 
research for at least two decades. 

B. THESIS OBJECTIVES 

This work should serve as the basis of what is a long 
term ongoing project at the Naval Postgraduate School. It is 
the theoretical and experimental studies of the scattering 
characteristics of targets and the development of radar 
target identification schemes. 

The main goal of this work is to calibrate the contin- 
uous wave (CW) scattering range. This is the first step 
toward studing target scattering characteristics for target 
identification. Calibration has been performed by using an 
assortment of perfectly conducting metallic spheres of known 
cross sections. The theoretical solution has been tested 
against experimental results for the spheres of different 
diameters and two orthogonal polarizations, horizontal and 
vertical. After calibration, the next step was a series of 
measurements on the right circular cylinders having very 
thin walls and finite lengths made of brass and stainless 
steel. They come in three diameters and several different 
length to diameter ratios. Measurements of magnitudes and 
phases of back scattering fields from the tubular cylinders 
were carried out and corresponding back scattering cross 
sections are obtained. Those results should serve as the 
basis for future works in the ongoing project. A right 
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circular cylinder is chosen because of its resemblance to a 
missile body. Progressively more complicated models evolving 
from a tubular cylinder to a missile will be constructed and 
studied. Each model which represents a step forward in 
increased complexity from a simpler one will be considered 
as a perturbation to the physical structure of the simpler 
model. By comparing the scattering data of the models to 
those of tubular cylinders, the effects of the successive 
perturbations to the physical structure on the back scat- 
tering cross section and phase shift can be observed. 
Explanations of these effects is currently being attempted. 
Application of this technique toward target identification 
and radar cross section reduction will be explored. 

In chapter II the general approach to the solution of 
the problem of electromagnetic scattering is given. In 
particular, the exact solution for the scattering by a 
perfectly conducting sphere is presented along with the 
considerations on the back scattering cross section. 

Chapter III describes experimental setup in the scat- 
tering laboratory of the Naval Postgraduate School and the 
measurement procedure, as well as the validation of measure- 
ment data. 

Computer programs for the evaluation of the data for a 
sphere, system calibration and target data acquisition are 
described in the chapter IV. 

Chapter V contains a summary and conclusions. It also 
discusses results and problems encountered during this work. 
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II. ELECTROMAGNETIC SCATTERING PROBLEMS AND SOLUTIONS 



A. GENERAL APPROACH TO THE PROBLEM 

A typical problem in electromagnetic scattering involves 
three principal elements; these are: 

1. A transmitting system consisting of an antenna or 
radiator that is excited by a generator of an electromag- 
netic disturbance of known behavior in time; 

2. An obstacle of arbitrary shape, size, orientation, 
and material that is fixed or is moving in any direction and 
manner at an unrestricted distance from the source; 

3. Receiving antenna with associated apparatus for 
determining amplitude, phase, and the polarization of the 
electromagnetic field at any point in space, including the 
surface (or surfaces) of the obstacle. 

The oscillating charges and currents in the driving or 
primary antenna constitute the source. The electromagnetic 
field maintainedby them is the primary or incident field. It 

~ ->i 

is denoted by a superscript "i": thus, E , H • The parasitic 
element under study is known as obstacle, radiator, scat- 
terer or target, depending partially on its shape and func- 
tion. The current in the source induces time-varying 
distributions of oscillating charges and currents in the 
scatterer. These, in turn, maintain an electromagnet i c field 
that is known as a scattered, or reradiated field. It is 
denoted by a superscript "s": thus, E > H • The total field 

-Hp -^p 

E » H , is the vector sum of the incident and the scattered 
fields. It is a resultant field at any point when the 
obstacle is present and the source is active. Thus, 

ftT £tT fti , £jS „ . r, 

E=E+E,H-H+H , Figure 2.1. 
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In general, the incident field is not independent of the 
scattered field in the sense that the currents and the 
charges in the source are not the same when the obstacle is 
present as when it is removed to infinity. The incident 
field does depend upon the presence and nature of the 



Z 




Figure 2.1 Configuration of a Scattering Problem 

obstacle, at an arbitrary point in space. In order to 
simplify the theoretical problem, it is usually assumed that 
the source is so loosely coupled to the obstacle, usually by 
a sufficiently large separation, that the currents and 
charges in the source may be regarded as independent of 
those in the obstacle. In this case, the incident field at 
any point may be treated as a quantity that is independent 
of the scattered field; it may, therefore, be determined in 
the absence of the obstacle. Alternatively, the scattered 
field may be obtained by subtracting the incident field in 
the absence of the obstacle from the total field at the same 
point when the obstacle is present. 
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The fundamental problem in scattering is the determina- 
tion of the scattered field amplitude, phase and polariza- 
tion. The scattering of EM wave from an obstacle is a 
boundary value problem. It is often hard, if not impos- 
sible, to obtain the exact solution, unless the shape of the 
obstacle coincides with one, of the limited, simple coordi- 
nate systems. 

To solve the problem, it is necessary, as usual, to 
start from Maxwell equations. The first step is to determine 
the scattered EM field in the surrounding space due to the 
currents on the surface of the scatterer. The boundary 
conditions on the surface are then applied to obtain a set 
of coupled equations, containing all components of the 
surface current. Only recently has become practical with 
the use of high-speed digital computers to estimate surface 
currents flowing on the body by a finite-difference solution 
to a network of simultaneous equations [Ref. 5]. Because of 
the limitations of computation time and storage capability, 
the finite-difference solution is applicable only when the 
dimensions of the target do not exceed a few wavelengths. 
For almost all practical problems, the use of approximate 
methods is the only practicable approach. 

With surface current distribution determined, other 
quantities such as the far field pattern and the scattering 
cross section are obtained from the scattered field evalu- 
ated at the begining. 

B. SCATTERING BY A SPHERE 

A sphere is the simplest object susceptible to exact 
solution by the method of separation of variables. The 
literature on electromagnetic scattering by spheres is 
extensive, and a great deal is now known in detail 
concerning the surface current distribution and the 
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scattered field in the neighborhood of the sphere [Ref. 6]. 
The exact solution of the scattering of a plane EM wave by 
an isotropic, homogeneous sphere was obtained by Mie in 1908 
and is usually refered as the Mie solution. 

1. Exact Solution 



The geometry of the problem in spherical coordinates 
r , q , $ is depicted in the Figure 2.2. 



Z 
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POINT 
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Figure 2.2 Geometry of Scattering by Sphere 



Let us consider a sphere with a relative dielectric 



constant e and an incident wave propagating in the 
z-direction and polarized in x-direction. 




( 2 . 1 ) 
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In the spherical coordinate system it is possible to 
express any electromagnetic field in terms of two scalar 
functions (Hertz potentials) which are the radial components 
of the electric and magnetic Hertz vectors [Ref. 7], 



n 

e 



III • 



r 



and 




( 2 . 2 ) 



Both IIi and II 2 satisfy a scalar wave equation 
2 2 

( V + k 2 ) • II = 0 for r > a (2.3) 



Electric and magnetic field can be expressed as 



) ^ /v 

E = VxVx(rllir) + ia>uo Vx (rH 2 r ) 

(2.4) 

^ A A 

H = -icucVx (rH 1 r ) + VxVx(rH 2 r) 



It can be shown that the general solution of the 
scalar wave equation in spherical coordinates may be 
obtained by a linear superposition of the particular solu- 
tions , thus 

00 n ( rr, ) 

rll = r • Z. Z IT J (2.5) 

n= 0 m= - n n 

= E E (c V (kr) + d x (kr ) } (P ^ (cos0 ) } 
n=o m=-n n n n n n 

(a cos(m^) + b sin(m$)} 
m m 
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The isotropic, homogeneous medium will be considered 
to be a dielectric so that the propagation constant k 2 is 
real, and the isotropic homogeneous sphere is characterized 
by a propagation constant kj which may not be complex. The 
ratio of these quantities defines the relative refractive 
index . 



ki mi k 0 mi 

m = = = 

k 2 m 2 k o m 2 



C 2 . 6 ) 



The incident field, given in the equation (2.1) can 
be derived from the two scalar functions IIi and II 2 and are 
expressed in the following spherical harmonic functions: 



rH 



i 

l 




” i n-1 C 2n+l) 
n_1 nCn+1) 



Y (k 2 r ) P 1 ( cos© ) cos$ 
n n 



mi 



- — ? (k z r)P l (cos0)sin$ 

ki n - L n(n+l) 



(2.7) 



In the same way, we write general expression for the scat- 
tered fields for r>a using arbitrary constants a and b : 

n n 



rill = — - — n f 1 — - 2n+ - 1} a C (k 2 r)p : (cos0)cos$ 

ki n(n+l) n n n 

( 2 . 8 ) 

rnf = — — — E.— ^ n+ - ^ b L (k 2 r)P 1 (cos0)sin$ 

ki n ~ 1 n(n+l) 
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It can be shown that by applying boundary conditions 

we can get four linear equations in the coefficients a , b , 

n n 

c , and d . These equations can be solved for that sets of 
n ’ n ^ 

coefficients. Only first two are of interest here, and 
these are given by 



a 

n 



T (a) •'?' (6) - m? (8) V (a) 
n n n n 

C Ca) V (B) - mY (B) •?' Ca) 

n n n n 



(2.9) 



b 

n 

where 



mT (a) y' (b) _ w (g)w' ( a ) 
n n n n 



ms (a) v' (B) - ? (B ) S ' (a) 
n n n n 



a = k 2 a = 2-ira/X = 2irm 2 a/X ( 



6 = kia = 2irmia/X 0 = ma 



The scattered fields E^ and E^ far from the sphere are given 
by 



E 0 ie 



• ik 2 r 



Si (© )sin$ 



k 2 r 



( 2 . 10 ) 



b 0 xe 



ik 2 r 



'0 



S 2 (0 ) cosO 



k 2 r 



where 
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{a tt (cos0) + b t(cos0) } 



n n 



n n 



( 2 . 11 ) 




(a x (cos0) + b tt (cos0)} 



n n 



n n 



P 1 (cos© ) 



— P^cosQ) (2.12) 

, « n 



TT (COS0) 



n 



T (COS0) 



n 



sin© 



n 



2 . Scattering Cross Section 

The radar cross section of a target is a quantita- 
tive measure of the ratio of power density in the vector 
signal scattered in the direction of the receiver to the 
power density of the radar wave incident upon the target 
[Ref. 8]. The vectoral nature of the scattering interaction 
requires specification of the transmitter polarization and 
receiver polarization with reference to target orientation 
in three dimensions. Transmitter operating frequency must 
also be specified. The radar or back scattering cross 
section is usually defined to be independent of range to the 
target, but it depends upon target shape and material, the 
angle at which the target is viewed, radar frequency, and 
the polarization of the transmitting and receiving antennas. 
Independence of the back scattering cross section upon the 
range to the target is valid under the far-field condition, 
i.e., when the target is sufficiently far from the transmit- 
ting antenna to justify the assumption that the incident 
wave is planar in the neighborhood of the receiving antenna. 
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The conventional definition of the scattering cross 
section in terms of the electric field can be, 



r s 2 

a(Q,$) = lim 4irR 2 |— — | C2.13) 

R-«° E 1 

where R = r is the distance from receiving antenna to 
target. Particularly for the back scattering $ = 0 , 0 = II , 
and for the perfectly conducting sphere |m|-*° • Hence, 



y’ (a) 
n 

£ ’ (a) 
n 



Ca) 

n 

C Ca) 

n 



Equation for S 2 C 0 ) yields, 



S 2 Cir) = -i p»-l_(n+l/2) 

n 1 C Ca)c' Ca) 

n n 

because of, 



(2.14) 



(2.15) 



Y (x) 
n 




x 



2 




5 (x) 
n 



CS)* 5 H (1 ,><x> 
l n+^ 



TT (-1) 
n 



C-l) n+1 



n(n+l) 

2 



x C-l) 

n 



, , >.n n(n + l) 
J 2 



Equation (2.10) now yields, 
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E s 

x 



5 



(2.16) 




ik 2 r 
= Eo — 

k 2 r 



Z , ( -1 ) n 



1 (n+1/2) 

? (a)? (a) 
n 



= 0 



Finally, back scattering cross section and the phase shift 
can be obtained as, 



a 



B 



2 r s , 2 

lim 4 ttR I — — : — | = 
R-h* E 1 



4lT 

k 2 



Z. (-l) n-1 (n+%) 

n= 1 

? (ka)? (ka) 
n n 



(2.17) 



where k 2 = k , a = ka 



L C-l) n-i (n+h) 

Phase = - arg{^ ; } (2.18) 

Z (ka)? (ka) 
n n 



For almost any finite target, an examination of the 
scattering cross section as a function of frequency leads to 
a separation of the frequencies into three regions each of 
which is characterized by a particular type of behavior 
[Ref. 9]. There is, first of all, the LF region where the 
cross section increases rapidly and smoothly with fourth 
power of frequency. This Reyleigh- region behavior is char- 
acteristic of any object whose greatest dimension is much 
smaller than a wavelength. Figure 2.3 shows the variation 
of radar cross section with the ratio of sphere circumfer- 
ence to wavelength, 2ira/X , where "a" is the sphere radius 
[Ref. 10]. 
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Increase in frequency, (or normalized circumference 

2ra/A, for sphere) introduces oscillations in the cross 

section curve which damp rapidly with increasing 2ira/A and 

2 

which vary about a mean value o/ia = 1 .The region where 
these oscillations are prominent C 1< 2ra/X< 10 ) is called Mie 
or resonant region, as mentioned earlier. Even for 
ka = 2ira/A = 43 cross section is still varying by ±1% about 




Figure 2.3 Radar Cross Section of the Sphere 
its mean value. 

In the high-frequency limit ( ka>10 ), known as 

optical region, the cross section settles down to the spec- 

2 

ular value o = ia , the projected area of the sphere. 
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III. MEASUREMENTS AND RESULTS 



The value of experimental radar cross section work is 
necessary and quite obvious when one realizes that we have 
exact solutions for the scattering from only a few bodies, 
although there are, of course, excellent approximate solu- 
tions for several simple shapes. Much emphasis is placed on 
approximate solutions in studying radar cross sections, and 
experimental results have been very important in helping to 
develop such solutions and in providing their validity. 

The principal difficulties experienced in measuring 
radar scattering cross sections are in the separation of the 
transmitted signal from the received signal and the elimina- 
tion of unwanted signals that are reflected from the back- 
ground, foreground, and the target model support. 

Radar cross section measurement facility at the Naval 
Postgraduate School utilizes CW system that depends entirely 
on cancellation the unwanted signals. The CW system, prob- 
ably the most reliable, is particularly useful for small 
models that are measured at ranges of 15 meters or less. 

This chapter deals with the description of the above 
system, its experimental setup, measurement procedures, as 
well as the validation of the obtained measurement data, in 
great detail. 

A. EXPERIMENTAL SETUP 

The configuration of the entire system is shown in 
Figure 3.1. System is divided into three basic parts as 
f o 1 lows : 
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1 . Range Geometries 



The physical arrangement of the transmitting and 
receiving antennas, target mount, walls, etc. is usually 




Figure 3.2 Anechoic Chamber at Naval Postgraduate School 

called the range geometry. A variety of range geometries 
have been developed: these are characterized by the way 
they have been designed to help eliminate unwanted signals 
reflected from the foreground and background [Ref. 11]. The 
range geometries commonly used are free field ranges with 
the antenna beams horizontal, oblique, or vertical; ground 
plane ranges, and image plane systems. Free space and image 
plane ranges are commonly used in both indoor and outside. 

The new "free-field" CW scattering range at the 
Naval Postgraduate School imploys a metallically shielded 
anechoic chamber Figure 3.2, [Ref. 14]. It can also serve 
as the free-field video impulse scattering range [Ref. 13]. 
Physical dimensions of the anechoic chamber are approxi- 
mately: length 6m, width 3m and height 3m. Anechoic chamber 
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has the very important advantage of being free from time 
varying interference due to external noise such as atmos- 
pheric noise, man-made noise (radio, TV, radar, etc.), as 
well as the weather conditions. This advantage is obtained 
at the expense of introducing critical problems that arise 
with reflections from the side walls, ceiling, and back 
wall. These reflections are reduced by lining the interior 
of the chamber with EM waves absorbers. EHP-8 microwave 
absorber made by Rantec Inc. is utilized which is designed 
for wide range of frequencies above 1.5 GHz. The units are 
fabricated of high quality urethane foam which is cut into a 
precise pyramidal configuration for front and back walls, 
and wedge configuration for lateral walls, floor and 
ceiling. Continuing improvement in absorber performance and 
chamber design and in electronic schemes for eliminating 
unwanted reflections have made it possible to make measure- 
ments on practically all types of targets except those that 
are too large to be compatible with the chamber dimensions. 



TABLE I 



Antenna Spe 

Frequency range 
Gain 

VSWR (maximum) 

Isolation 



if icat ions 

4-18 GHz 
6 - 12 dB 
3.2 : 1 

<20 dB below 5.5 GHz 
>20 db above 5.5 GHz 



The target is supported by a low cross section Styrofoam 
column to minimize reflections [Ref. 12]. Detailed analysis 
of the support influence is given in section "C" of this 
chapter . 

The Naval Postgraduate School radar cross section 
range operates as bistatic system utilizing two horn 



28 



antennas. Both transmitting and receiving antennas are 
mounted on a removable panel in the front wall. Base 
aluminum plate that bears both antennas has mechanism 
allowing antenna adjustment in all three axes. Adjustment 
screws are accessible from outside the anechoic chamber 
providing beam steearing toward target. Antennas are 
designed for both, linear polarizations, horizontal and 
vertical. While used as horizontally polarized, transmitter 
and receiver, the other feeds are terminated with matched 
loads and vice versa for the vertically polarized feeds. 
Basic antenna characteristics are given in the Table I. 

2 . Instrumentation 

Equipment used in the scattering laboratory of the 
Naval Postgraduate School, shown in the block diagram Figure 
3.1, are listed below: 



• 


Microcomputer 


HP 


-85 




• 


Synthesized Signal Generator 


HP 


8672A 




• 


Harmonic Frequency Converter 


HP 


8411A 




• 


Network Analyzer 


HP 


8410C 




• 


Phase -Magnitude Display 


HP 


8412B 




• 


RF Amp 1 i f i e r 








• 


DC Power Supply 


HP 


6227B 




• 


Digital Multimeter 


HP 


3466A 




• 


Digital Voltmeter 


HP 


3456A 




• 


Digital Voltmeter 


HP 


3455A 




• 


Flexible Disc Drive 


HP 


82901M 




• 


Plotter 


HP 


7225B 




• 


HP Interface Bus 








• 


Directional Coupler 








• 


Coaxial Cables 










HP-85 microcomputer provides 


necessary 


signal 


processing using developed computer codes, 


controls 


Hp 8672A 
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signal generator, both HP 3456A and HP 3455A digital voltme- 
ters, as well as the HP 82901M disc drive and HP 7225B 
plotter. Communication between microcomputer and remote 
controled devices takes place through HP Interface Bus 
(HP- IB ) . 

The HP Model 8672A synthesized signal generator 
produces RF oscillation in a frequency range of 2 to 18 GHz. 
The output is leveled and calibrated from +3 to -120 dBm. AM 
and/or FM modulation modes, and most other modes can be 
remotely controled using the HP-IB programming format. 

Rf signal from signal generator enters the microwave 
amplifier whose gain is 23 dB . Because of the low level 
output from signal generator, amplification is necessary in 
order to obtain satisfactory S/N ratio in the receiver 
output of the target return at higher frequencies. The power 
supply (HP 6227B) provides the amplifier with a DC voltage 
and the digital multimeter (HP 3466A) monitors that voltage. 

Passing through directional coupler, the Rf signal 
is split into two. Unattenuated RF portion feeds transmit- 
ting horn antenna and the other portion is attenuated by 43 
dB and feeds the reference port of the harmonic frequency 
converter. The test port input of the converter is fed by 
the return signal collected with the receiving horn antenna. 

The HP 8411A harmonic frequency converter and the HP 
8410C network analyzer convert the RF signal being measured 
to 278 kHz signals that have the same magnitude and phase 
relationships. These two units together with the HP 8412B 
phase-magnitude display, function as a phasemeter and a 
ratiometer for direct, continuous, simultaneous phase and 
magnitude ratio measurement of RF voltages. The complete 
network analyzer measures phase angles 0 to 360 and magni- 
tude ratios in decibels over a dynamic range of 60 dB . 
Harmonic frequency converter automatically tracks the 
frequency of the signal applied to the reference input. 
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The magnitude and phase DC voltages, taken from the 
phase-magnitude display unit, are measured by the HP 3456A 
and HP 3455A digital voltmeters respectively. Those are 
versatile digital voltmeters with AC, DC, OHMS, and math 
functions, and offer the capability of remote programming. 

•3 . Targets 

Targets used in this measurement are several 
aluminum and steel spheres of different diameters and series 
of right circular cylinders of various lengths and diame- 
ters, made of thin walled brass and stainless steel. Spheres 
are used for calibration as canonical or standard targets of 
known (computed) radar cross sections. This measurement 
reveals correlations of the theoretically obtained results 
using the Mie series given by the equations (2.17, 2.18). 

Cylinders, as primitive models of the missiles, are used to 
form necessary data base for future work. Adding caps, fins 
and wings, this primitive model will be perturbed toward a 
real missile body. Exhaustive measurements of radar cross 
section and phase shift will be taken out, and results 
should reveal versatility of this scheme in its application 
to target identification and radar cross section reduction. 

One advantage of the radar range, as mentioned 

earlier, is the practicability of testing models that are 

smaller and cheaper than full-scale targets. For this 

purpose the radar wavelength is scaled by the same factor as 

the dimensions of the model. If "D 0 " is any given dimension 

of the target, and "D " is the equivalent dimension of the 

M 

model, the following scaling relation is employed, 




C 3 . 1 ) 
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where X is the wavelength used for the measurement and X . 
is the wavelength for which the target cross section is 
required. At the same time, the measured cross section of 
the model should be multiplied by p 2 to give the target 
cross section a Q , 



a 



o 




p 2 a 



M 



(3.2) 



Exact scaling requires the model conductivity to be equal 
the target conductivity multiplied by the ratio ”p" . This 
parameter is not critical if the model conductivity is high 
enough . 



B . MEASUREMENT PROCEDURE 

Procedure for measuring back scattered electric fields 
from a target was suggested in chapter II. Metallic spheres 
are commonly used as standards since their cross section can 
be computed, equation (2.17). Prior to the start of taking 
measurements on various targets, system transfer function 
has to be obtained through a calibration program. 

1. Calibration of the System 

After run of the calibration program, microcomputer 
HP-85 remotes and initializes all programable units, Figure 
3.1. Signal generator is set on and begins to generate RF 
energy. This RF , CW signal undergoes 23 dB amplification in 
order to improve S/N ratio. It passes directional coupler 
and is radiated, by the transmitting antenna, toward the 
target into the anechoic chamber. Signal return at the cms 
receiving antenna, which is the vectoral sum of the target 
echo and the background radiation is picked up by the 
receiving antenna which feeds the "test" port of the 
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harmonic frequency converter. The background radiation 
includes the direct coupling between the two antennas and 
the back scattered signal from the chamber. The "reference" 
port of the same instrument is feed with the RF signal 
passed through another branch of the directional coupler. 
This signal is attenuated by 43 dB and enters the frequency 
converter at the about -20dBm. The converter translates the 
RF signals being measured to 278 kHz while preserving the 
same magnitude and phase relationships. It also tracks auto- 
matically any change in frequency applied to its reference 
port, and passes that signal to the network analyzer. The 
frequency converter, the network analyzer together with the 
phase-magnitude display unit function as a phasemeter and a 
ratiometer for direct, continuous, simultaneous phase and 
magnitude ratio measurement of RF voltages, as mentioned 
earlier. Both magnitude and phase of the signal return are 
measured by the HP 3456A and HP 3455A digital voltmeters 
respectively, and stored in the HP-85 microcomputer which 
then computes the real and imaginary parts of the signal. 
Separation of the target echo from the background radiation 
can be accomplished by measuring the background radiation 
data as will be explained latter. 

The procedure repeats "NO” times at each frequency, 
and after the average value of magnitude and phase are 
computed, the signal generator increments frequency by one 
step, and the system goes through new cycle. The block 
diagram of Figure 3.3 presents various transfer functions of 
the scattering range [Ref. 13], [Ref. 14]. Those are: 

• H t ( to ) - Transfer Function of the Transmitting Circuit 

• H (oi ) - Transfer Function of the Receiving Circuit 

K. 

• h a<“ ) - Antenna Coupling 

• H c (to ) - Clutter Signal From the Anechoic Chamber 
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• Hg(oj) - Back Scattered Signal Normalized ( |H ( w ) | 2 is 

the cross section and arg[H (co)] is the phase shift of 

the target) 

T 

• H g (w) - Back Scattered Signal of an Arbitrary Target 

Besides, X(oa) is input signal provided by the signal gener- 
ator and the RF amplifier. Y(oj ) is reading of the network 
analyzer. The network analyzer provides the ratio Y(co) 'of 
the test channel signal to the reference signal which is a 
complex number at each frequency step. Note that every 
block, including X(w) contains noise. 



X(w) 



Y 3 Coo) 




Figure 3.3 System Representation with Target in Place 



back 

rate 



The system transfer function H (oj) and the target 
T sys 

scattering H g (co ) can be found by performing three sepa- 
measurements : 

1. without the target in place 

2. with the canonical target in place 

3. with the arbitrary target in place 



From the first measurement, background data (response from 
the chamber without target in place) are obtained. Block 
diagram of Figure 3.4 represents this measurement. 
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XCoo) 



Yj C co 




Figure 3.4 System Representation in Background Measurement 

Signal return Yi (u) which corresponds to the scat- 
tered signal of the background can be expressed by the equa- 
tion (3.3). 



Yi(u) = H (co) { H (co) + H ( co ) }H (oo) C3.3) 

T A C R 



The procedure described above is utilized in the 
second measurement as well, in order to acquire data with 
the sphere as canonical target in place. For this case 
system response is labeled as Y 2 ( co ) and the measurement is 
represented with the block diagram of Figure 3.5. Equation 
(3.4) represents this configuration. 



These two measurements provide system calibration 
because system transfer function can be found: 
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(3.4) 



Y 2 C Cjo ) = H t (w) • {H a C cu ) + HgCu) + H c Cco)} -H (a>) 



Subtracting background given by equation (3.3) from the data 
obtained from canonical target given by equation (3.4) 



XU) 



Y 2 Coj 




Figure 3.5 System Representation with Canonical Target 

yields scattered data from the sphere (reference standard) . 
Let it be Y 4 ( oj ) : 



Y 4 Co)) = Y 2 U) - YiCoj) = H CoOHgCciOH Cw) 



( 3 . S ) 



Q 

In equation (3.5) H (co) can be 

S 

equation (2.17). Equation (3.5) 
function of frequency "oi" : 



obtained from Mie solution, 

yields H (co) as complex 
sy s 
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-1 



C 3 . 6 ) 



H (co) = { H (a) ) H (co) } 
sys T R 



= (co ) Yu (co ) _1 



After the calibration of the system is completed and 

its transfer function H (to) obtained, we can find the 

sys 

scattered field for any arbitrary target. To accomplish 
this, we have to perform a third measurement with the target 
in place. System representation is given on the block 
diagram of Figure 3.3 and its mathematical interpretation is 
as follows: 



Y 3 Cco) = H (to) (H (co) + H (co ) + K^(co))H (gj) 

I ii C SR 



(3.7) 



The background data Y x (co) are subtracted from the target 
data Y 3 (co) which we call Y c (co), 



Y s Cco) = Y 3 (co ) - YiCco) = H Coi)H^Caj)H (to) 



(3.8) 



T 

Target transfer function H (co)is obtained by simply pluging 

s 

equation (3.6) into equation (3.8). It yields : 



H^(co) = Ys (co) H (oj) 
* sys 



C3 . 9 ) 



Radar or back scattering cross section can be computed as : 
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(3.10) 




2 



where | ( o> ) | 2 = o 

& 



TH 



is theoretically computed back scat- 



tering cross section of the sphere used as reference stan- 
dard. The phase shift of the target is : 



2. Measurement on Spheres and Cylinders 

Measurement procedures for both calibration and a 
target data acquisition are similar as was described in 
previous section and consists of three steps. After calibra- 
tion data from several spheres are obtained and compared 
with the computed results of the respective target. For that 
purpose, theoretical results were obtained using Mie solu- 
tion, given by the equation (2.17), and represented in 
graphical form , Figures 3.7 to 3.13. Computer plots of 
target cross section and phase were taken for three spheres 
of different diameters, in the frequency range of 2 to 18 
GHz, incrementing frequency in steps of .02 GHz. For the 
first part of the experimental work, the targets used were 
spheres and in the second part were cylinders. All targets 
are labeled with numbers 1 through 8 for convenience. 
Besides, brief description of each target is given in the 
header of the corresponding graph. It contains basic phys- 
ical dimensions of the target, and corresponding measurement 
specifications such as the polarization of the antennas and 
the sphere used for calibration. A brief review of the 
targets characteristics is given in Table II. All graphs of 




(3.11) 
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theoretical computations and experimental results are given 
at the end of this chapter. 



TABLE II 

Target Characteristics 



TARGET 1 


Sphere , 


made 


of 


alum. , 


D 


= 


.875" . 




TARGET2 


Sphere , 


made 


of 


steel , 


D 


- 


1.125" . 




TARGET3 


Sphere , 


made 


of 


steel , 


D 


= 


3.187". 




TARGET4 


Cylinder , 


made 


of 


brass , 


D 


= 


•75" , 


L = 3.875". 


TARGET5 


Cylinder , 


made 


of 


steel , 


D 


= 


•75" , 


L = 3.875". 


TARGET6 


Cylinder , 


made 


of 


brass , 


D 


= 


• 5" , 


L = 3". 


TARGET 7 


Cylinder , 


made 


of 


brass , 


D 


- 


.375" , 


L = 1.875". 


TARGET8 


Cylinder , 


made 


of 


brass , 


D 


= 


• 75" , 


L = 6". 



Measurements on spheres and cylinders were done for 
a frequency range of 10 to 15 GHz in steps of 0.1 GHz. 
Graphs of target cross section and phase versus frequency 
are represented by dots overlayed with corresponding theo- 
retical curves. 

C. VALIDATION OF MEASUREMENT DATA 

Recognition of sources of errors in measurement data is 
usually difficult. In a few special cases, it is possible to 
recognize the presence of error and determine its source by 
observing the deviation of a target cross section vs. 
frequency plot from anticipated behavior. The most common 
difficulties are system noise and target-range problems. The 
latter include 1). the requirement that the field incident 
on the target (model) be a reasonable approximation to a 
uniform plane wave over the volume of space occupied by the 
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scattering object, 2). the interference of extraneous 
reflections from objects other than target, and 3). method 
of supporting the target [Ref. 12], [Ref. 15]. Each of 
these problems as well as their effect to the measurement 
results will be discussed separately. 

1 . System Noise 

The main sources of the system noise come from the 
RF amplifier, the network analyzer and from thermal noise in 
the Rf signal path and the electronic circuitry. By using 
the anechoic chamber which is metallically shielded, almost 
all external noises were eliminated. The noise from the Rf 
amplifier is present in Y(co) because the test channel signal 
is a combination of signals output from the amplifier at 
different instances while the reference channel signal left 
the amplifier at a particular instance. Thus their ratio 
will contain time varying noises. The network analyzer is a 
noise source because for every set of target data obtained, 
the generator steps through the frequency several times. At 
one frequency, the frequency converter may not pick the same 
harmonics each time the generator steps on this 
frequency . This causes variations in the network analyzer 
responses. The other noises are thermal noises. All these 
components of the system noise sources produce random noises 
and can be significantly reduced by averaging. Averaging 
process was done twice, first during the acquisition of the 
magnitude and phase data at a single frequency. The entire 
measurement procedure is then repeated three or four times 
and averaged . 

2 . Antenna -Target Distance 

In general, the fields scattered by an object depend 
on the variation, in phase and magnitude, of the incident 
field over the surface of the object as well as on the 
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object's material, size, and shape [Ref. 16]. The scat- 
tering cross section of a target is defined for a plane 
incident wave and for the far field scattered by the target. 
Because of the confines of an indoor range, non-plane wave 
illumination of the target is a significant source of meas- 
urement error. The far from an antenna is usually spherical 
along its boresight axis and can be only an approximation to 
a plane wave. Such a spherical wavefront, when referenced to 
a plane perpendicular to the antenna boresight and tangent 
to the wavefront at the nominal target position, will 
produce, to the lowest order of approximation, a phase devi- 
ating quadratically from that of an ideal plane wave along 
directions perpendicular to the central ray. The plane wave 
approximation to the far field is satisfactorily estab- 
lished, as the rule of thumb, for ant enna- t arget distance 
equal to or greater than 2D 2 /A, where D is the largest 
dimension of the target (model). This condition as usually 
imposed relates to the uniformity of the phase front over 
the lateral region occupied by the target. To preserve both 
magnitude and phase uniformity within the test region, the 
magnitude gradient along the antenna boresight should be 
reduced, as much as possible. Due to reciprocity, the same 
rule of thumb of 2D 2 /A as the range for far field applies to 
the scattered field from the target . To avoid these diffi- 
culties, a more stringent far-field condition should be 
applied. In other words, target should be moved farther from 
the transmitting antenna. 

Initially measurements were done with antenna to 
target range of 1.5m, that caused disagreement with theory 
for the spheres 6" in diameter and larger. Distance antenna- 
to-target was changed to 3m, but the influence of the back- 
ground, that will be explained latter, as well as S/N ratio 
caused erroneous results. Finally, an optimal range of 2m 
was established to satisfy most of, often contradictory 
requirement s . 



41 



3 . Background Radiation 



Accurate target cross section measurements require 
the background signal to be held at very low level. 
Background radiation can be defined as returns from objects 
in the vicinity of the target. Background effects may be 
classified as direct backscatter and target scatter. Both 
sources of error are reduced by arranging radar- absorbing 
material so as to minimize the undesired background signal 
intensity. Residual direct backscatter can be determined, 
as it was done in this experiment, by measuring the return 
in absence of the target . Target scatter, or bistatic 
coupling effects frequently arise from the strong target- 
scatter lobes in the forward direction (away from the 
radar) . They characterize many targets and are often two or 
three orders of magnitude greater than the direct-back scat- 
tering which we measure. Tor instance, a sphere which diam- 
eter equals 6” has aforward- scatter about 25 dB stronger 
than the back scattering cross section [Ref. 17]. Thus the 
possibility that the scattered energy from the target may 
illuminate the surroundings to return significant energy to 
the radar cannot be ignored. Illumination of the background 
by target- scattered energy may be followed by reradiation 
direct to the radar or by scattering back to the target and 
then to the radar. 

This target- scatterer effects were not measured in 
the experiment. Even direct backscatter can not be deter- 
mined quite accurately because the introduction of the 
target alters the radar illumination of the background as 
well, so the background cancellation is no larger complete. 
The resulting error depends upon the target cross section, 
and the magnitude of the background before and after the 
target introduction. It is generally found that the error is 
small if the radar- to - target distance is small compared with 



42 




Figure 3 . 6 Measurement Error Due to Background 

both the radar- to-background and target- to-background 
distances but that it can become large when the target 
distance is increased. As mentioned earlier, background 
errors were noticed by taking measurements at distances of 
1 . 5m and 3m . 

Figure 3.6 shows the measurement error that can 
result from the background as a function of the ratio of 
effective background cross section to target cross section 
[Ref. 18]. To prevent error greater than ldB , that corre- 
sponds to 20% in cross section, it is necessary that the 
effective background cross section be less than one 
hundredth the cross section of the target. 
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4 . Target Supports 



Target support, as a structure employed to control 
and permit measurement of the target position and orienta- 
tion, is another important source of measurement error. The 
magnitude of the error introduced depends upon the back 
scattering from the target support relative to that of the 
target. The error is significant when the target cross 
section is small. If interference with the support is 
suspected, reruns of the experiment with different support 
configurations frequently reveal changes in both cross 
section and phase. In the case of the measurements 
performed on spheres, when the exact solutions and corre- 
sponding graphs for the target cross section and phase are 
known, deviations from those known patterns pinpoint to the 
errors in the measurements. 

Second effect that may occur is modification of the 
incident field by the support device. This effect, since it 
represents coupling between the target and its support, is 
impossible to evaluate quantitatively. 

In this work, the effect of the back scattering from 
the support, as well as the effect of the field modifica- 
tion, were significantly reduced using tapered Styrofoam 
supports . 

Finally, a brief comment on measurement results will 
be given, along with computer plots of target cross section 
and phase for all targets, spheres and cylinders. 

Measurements on spheres were done to obtain system 
calibration and reveal its validity for three spheres of 
0.875", 1.125" and 3.187". Besides these measurements have 
shown almost all previously discussed effects on the accu- 
racy of the measurements itself. 

Figures 3.13, 3.14 shows target cross section and 
phase respectively of the 3.187" sphere in diameter. 
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Immediately after calibration, that same sphere was used as 
the target. Agreement between measurement and theory is very 
good, remembering that the effect of the system noise 
couldn't be eliminated completely as well as the target- 
scatter background. It is observed that the results change 
in time because of the random nature of the system noise, 
and thence the background changes. In the Figures 3.15 to 

3.18 same calibration was used and results are not as good 
as in previous measurement, but they still are within ±8% of 
the theory. The smaller the target cross section the bigger 
the deviation. This is in accord with the consideration of 
the background, S/N ratio and support influence on measure- 
ment results . 

When the calibration was performed using 1.125" 
sphere, and the data taken of that sphere and 0.875” sphere 
respectively, disagreement between measurement and theory 
became more pronounced. In this case, smaller cross section 
of the target decreases S/N ratio, background and support 
effects are more enhanced, that all contributes to the 
bigger measurement error. Results are presented in Figures 

3.19 to 3.22. 

Influence of the target support reflections and 
coupling are revealed in Figures 3.23, 3.24 in comparison to 
the Figures 3.13, 3.14. This time TARGET3 (3.187" sphere) 
was used for calibration and to obtain target data as well, 
but with the cylindrical support instead of original one 
that was tapered. Both cross section and phase are changed 
slightly. If the target was sphere of smaller cross section, 
this effect would have been more noticable. 

Measurements on cylinders were done using 3.187” 
sphere for system calibration (CAL3). Same frequency range , 
as for the spheres, from 10 to 15GHz was chosen, with 
frequency increments of 0.1 GHz. Measurement results are 
plotted with solid lines to reveal approximate shape of the 
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target cross section and phase shift. Results obtained in 
several different measurements, separated in time, were 
reproducible and thence could be basis for future work on 
the perturbation of a cylinder to a missile body. 

Figures 3.27 to 3.30 show radar cross section and 
phase for two cylinders of the same dimensions, but made of 
two different materials, brass and stainless steel. 
Objective was to show influence of conductivity on the 
changes in target cross section and phase shift. This effect 
couldn't be observed because the both targets are made of 
materials with comparable , very high conductivity. 

In the Figures 3.31 to 3.36 targets were cylinders 
of different diameters and lengths. Corresponding graphs 
show that the cross section patterns depend upon cylinder 
diameter only (number of resonances is the same), while the 
target length along with its diameter contributes to the 
absolute value of the cross section. Phase on the other hand 
depends only upon the cylinder diameter. 

Changed polarization from horizontal to vertical 
gave quite different results for TARGET4 and TARGET6 in 
comparison to those of horizontal polarization obtained for 
the same targets, as was expected. Figures 3.36 to 3.40 
present these results. 
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THEORY. SPHERE □=. 875 inch 
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Figure 3.7 TARGETl - Cross Section (Theory) 
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FREQUENCY CGHz) 

Figure 3.8 TARGET1 - Phase (Theory) 
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FREQUENCY CGHz) 

Figure 3.9 TARGET2 - Cross Section (Theory) 
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FREQUENCY CGHzD 

Figure 3.10 TARGET2 - Phase (Theory) 
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FREQUENCY (GHz) 

Figure 3.11 TARGET3 - Cross Section (Theory) 
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FREQUENCY (GHz) 

Figure 3.12 TARGET3 - Phase (Theory) 
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FREQUENCY CGHz) 

Figure 3.13 TARGET3 - Cross Section 
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Figure 3.14 TARGET3 - Phase 
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Figure 3.15 TARGET 1 - Cross Section 
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Figure 3.17 TARGET2 - Cross Section 
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Figure 3.19 TARGET2 - Cross Section 
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Figure 3.22 TARGET1 - Phase 
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Figure 3.24 TARGET3 - Phase 
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FREQUENCY (GHz) 

Figure 3.25 TARGET3 - Cross Section 
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Figure 3.26 TARGET3 - Phase 




<Sd3JL3W ’OS) 



NOI133S SSOdD 



67 



Figure 3.27 TARGET4 - Cross Section 
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Figure 3.28 TARGET4 - Phase 
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Figure 3.30 TARGET5 - Phase 
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Figure 3.31 TARGET6 - Cross Section 
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FREQUENCY (GHz) 
Figure 3.32 TARGET6 - Phase 
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Figure 3.34 TARGET7 - Phase 
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Figure 3.38 TARGET4 - Phase 
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Figure 3.40 TARGET6 - Phase 



IV. DESCRIPTION OF COMPUTER PROGRAMS 



The computer programs described in this chapter are 
written in BASIC, for frequency domain scattering measure- 
ments at the Naval Postgraduate School. They are run on the 
Hewlet Packard - 85 microcomputer which, along with HP 
82901M Flexible Disc Drive and HP 7225B Plotter, controls 
programable instruments in experimental setup (Figure 3.1), 
performs necessary computations and store and plot results 
of the measurements. Four programs were written to complete 
scattering experiment. First, "SPHERE . DRIVEO" program 
computes theoretical values for perfectly conducting sphere 
which are latter used for system calibration by the program 
"CALIB . DRIVEO" . When the calibration is completed and 
system transfer function obtained, program "TARGET . DRIVEO" 
serves for target data acquisition and computation of its 
back scattering cross section and phase shift, and to plot 
those results immediately after measurement completion. 
Measurement results can be stored in the files on floppy 
disc, and if we want them to plot without taking measurement 
again, it is convenient to use separate program 
"PLOT . DRIVEO" . Programs are written to be interactive and 
user friendly. All that is necessary to run any of those 
programs is to load them into main memory and then issue the 
command "RUN" . 

A. THEORY OF SPHERE 

This program, under the name " SPHERE . DRIVEO" , computes 
theoretical values of the back scattering cross section and 
the phase shift of a perfectly conducting sphere, using Mie 
solution given by equations (2.17, 2.18). Theoretical values 
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are stored in files of 800 records, one for each frequency, 
in the frequency range from 2.02 GHz to 18 GHz, at 0.02 GHz 
steps. The main part of the program (lines 10 through 660) 
creates the file named "H$" , that stores one value of the 
cross section and one value of the phase, for each frequency 
increment of 0.02 GHz. It calls in SUBROUTINE I (lines 710 
through 990), and this subroutine by itself calls in 
SUBROUTINE II (lines 1040 through 1530), which together 
carry out the task of computation using previously mentioned 
equations (2.17, 2.18). Program is listed in following 
pages . 
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10 
20 
30 
40 
50 
60 
70 
80 
90 
1 00 
1 10 
1 20 
130 
140 
150 
1 60 
170 
18© 

1 90 
200 
210 
220 
230 
240 
250 
260 
270 
230 
290 
300 

310 HT=" THEOR 1 . DRIVE 1 " 

320 ! 

330 A0=1 . 125*. ©254/2 ! SPHERE 
RADIUS IN METERS 

340 ! 

350 ! 

360 X9=2*PI ! PARAMETER 

370 ! 

380 ! 

390 01=2 ! STARTING FREQ IN GHZ 
408 02=18 ! FINAL FREQ IN GHZ 
410 Q4= . 02 ! FREQ STEP IN GHZ 
420 ON ERROR GOTO 440 
430 PURGE H $ 

440 OFF ERROR 

450 CREATE H*, 800, 16 ! OPEN A 
NEW FILE WITH 80Q RECORDS 
460 ! OF 16 BYTES EACH. EVERY 
470 ! RECORD STRORES ONE MAGNI 
TUDE AND ONE PHASE DATA. 

480 ! 

490 ASSIGN# 1 TO Hf 
500 DIM BSC 1 44 > , B9 C 1 44 > , D3 C 1 44 > , 
D9 <144) ! 1 44>L8 = INT <2*K0*A0 

+ 3 > 

510 F0=Q 1 

520 FOR 1=1 TO 800 



' "SPHERE . DRIVE© 



PROGRAM TO COMPUTE 
VALUES OF BACKSCRTTER I NG 
CROSS SECTION AND PHASE 
FOR A PERFECTLY 
CONDUCTING SPHERE. 



THEORETICAL VALUES ARE 
STORED IN FILES OF 800 
RECORDS, ONE FOR EACH 
FREQUENCY FROM 2.02 GHZ 
TO 18 GHZ hT .02 GHZ 
STEPS . 



"THEOR© . DRIVE 1 " FOR THE 
. 875 11 DIAMETER SPHERE 
"THEOR! . DRIVE 1" FOR THE 
1.125" DIAMETER SPHERE 
"THE0R3 DRIVE1" FOR THE 
3.187" DIAMETER SPHERE 



FILE HT STORES THE 
COMPUTED RESULT 
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530 DISP "FREQ LOOP=",I 
S40 F h = F 0 + Q 4 

550 DISP "FREQ <GHZ>=" .■ F0 

566 Kl=.3/F0 ! WAVELENGTH 

570 K0=X9'K 1 ! WAVE HUMBER 

580 GOSUB 700 

59© DISP ”E E0 

600 DISP “P ="/P0 

610 PRINT# 1,1 ; E0,F'0 

620 NEXT I 

63© ASSIGN# 1 TO * 

640 CLEAR 

650 DISP "END OF COMPUTATION" 

660 END 
670 ! 

680 ! 

690 ! 

700 ! 

710 ! SUBROUTINE I 
720 1 

730 L0=INT<2*K0*A0+3> 

740 IF L 0 < 1 4 5 THEN 766 
750 DISP "K0*AO TOO LARGE FOR 
CURRENT ARRAY DIM" 

760 Z=K0*A0 
770 GOSUB 1000 
780 E8=0 
790 E9=0 

800 FOR N= 1 TO L0 
810 L=L0-N+ 1 
820 M8=D8 C L > •'"•2 + 09 < L '> A 2 
830 M9=B 8 < L > " 2-*-B9 < L > •'-2 
840 R7= < L+ . 5 > ''M8 -' M9 
850 A8=A7* < B9 < L > *D9 < L > -B8 < L >*D8 < 
L > ) 

860 89=07* < B8 < L > *09 <L>+B9<L> *08 < 
L > > 

870 E8=A8-E8 

880 E9=A9-E9 

890 NEXT N 

900 E 8 = - E 8 

910 E9=-E9 

920 E0=E8 x -2+E9'-2 

930 P0=ATN2<E9, E8) 

940 E0=E0/K0*2*K1 i CROSS-SECT. 

950 P0=P0-X9* I NT < P0/X9 ) 

960 IF P0 <P I THEN 980 
970 P0=P0-X9 

980 p0=-p0 

990 RETURN 
1 000 ! 

1010 ! 

1 020 ! 

103R 1 

1040 ! SUBROUTINE II 
1050 ! 
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1060 IF 2>L0- 1 THEN 1370 
1070 Z 2 = Z 2 2 
10:50 N2 = 2TZ2+L0+ 1 
1090 Dl=2*N2+3 
1100 02=01* <2*N2+5> 

1110 D3=D2*<2*N2+7> 

1120 D4=D3*<2*H2+9> 

1 130 F 1 = 1 -Z2^D1 +Z2' N 2--'< 21:02?' -22^3 
✓<6*D3> 

1 140 F2=2t <. 1 -'D 1 -22 ''D2+Z2^2^ < 2£03 
>-22 3 s < 6 1 D 4 > > 

1 150 N=2*Z2 
1166 S 1 =F 1 

1170 Fl=<2*M+l>*Fl/2-F2 
1180 F2=S 1 

1190 IF ft E S ■' F 1 > < 1 . E 1 0 0 THEN 1230 

1200 Fl=Fl*i . E-108 

1210 F 2 = F 2 T 1 . E-100 

1220 SI =S 1 T1 . E-100 

1230 M=M- 1 

1240 IF M + 1 > L 0 THEN 1160 

1250 B8<L0>=F2 

1260 B8<L0-1 >=F 1 

1270 N0=L©-2 

1280 FOR K= 1 TO N0 

1290 N=L0-K-1 

1300 B8<N>=<2*N+3>*B8CN+1 >^Z-B8< 
N + 2 > 

1310 NEXT K 

1320 R1 = <SIN<Z>/'Z-C0S<Z> B8k 1 > 
1330 FOR K= 1 TO L0 
1340 BS<K>=fll*B8<K> 

1350 NEXT K 

1360 GOTO 1420 

1370 B 8 < 1 > = S I N C 2 > / 2 - C 0 S < Z > 

1 380 B8 < 2 > = < 3-'Z~2- 1 > *S I N < 2 > -3*C0 
S < 2 > -'2 

1390 FOR N=3 TO L8 
14O0 B8 < H> = < 2TN- 1 > TBS C N- 1 > •'2- BSC 
N-2 ) 

1410 NEXT N 

1420 B9 < 1 > =-S I N C 2 > -COS < 2 > .'2 
1430 B9<2> = < 1-3'Z'-2>TC0S<Z>-3*SI 
N<2>/2 

1448 FOR N=3 TO L0 
1450 B9 C N> = <2TN- 1 >*B9<N-1 >^Z-B9( 
N-2 > 

1460 NEXT N 

1470 D8< 1 > = <1-1/'Z A 2>*SIN<Z>+C0S< 

2 > / Z 

1480 D9<T> = < 1, '2^2-1 >*C0S<2 > + S IN ( 
2 > ,'Z 

1490 FOR N=2 TO L0 
1 500 08 < N > =B8 •: N - 1 > -N*B8 ( N > ,'2 
1510 D9(N> = B9(N-1 >-NTB9<N>.'Z 
1520 NEXT N 
1530 RETURN 
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B. CALIBRATION PROGRAM 



The calibration program "CALIB . DRIVEO" computes system 
transfer function using a perfectly conducting sphere 
(canonical target), over the frequency range of 10-15 GHz, 
at 0.1 GHz steps. Computation is based on theoretical values 
using the program "SPHERE . DRIVE" , and the results are stored 
in the file "C$" . 

The main program (lines 10 through 1240) starts with 
comments and necessary declarations of the variables used 
(lines 10 through 520). Further, it takes theoretical 
values from H$ file and stores them into main memory. For 
the documentation of the results and user convenience, 
HEADER SUBROUTINE is called in (lines 2160 through 2280). 
To compute system transfer function, background data as well 
as the data acquired from canonical target must be avail- 
able. The main program ofers two possibilities, whether to 
use the background data obtained in the most recent calibra- 
tion, or to take them immediately after program is run. If 
user decides to take new background data, program remotes 
all devices establishing their initial setup (lines 530 
through 870), and calls in BACKGROUND DATA COLLECTION 
SUBROUTINE. This subroutine serves to increment frequency of 
the signal generator (device 719) Ml=51 times in steps of F9 
=0.1 GHz, and in each loop to take averaged data of the 
magnitude and phase, acquired by the digital voltmeters. It 
computes their real and imaginary parts, to be stored in the 
file A$ as well (lines 1300 through 1490). To do so, it 
calls in SUBROUTINE TO ENTER MAGNITUDE AND PHASE DATA from 
digital voltmeters. HP 3456A digital voltmeter (device 720) 
reads magnitude, and HP 3455A digital voltmeter (device 722) 
reads phase (Figure 3.1). In each loop both digital voltm- 
eters simultaneously take N0=10 readings of the magnitude 
and phase, average them and pass those values to the back- 
ground data collection subroutine (lines 1550 through 1830). 
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The main program proceedes to compute canonical target 
data, after it is placed into anechoic chamber, by calling 
in TARGET DATA COLLECTION SUBROUTINE . This subroutine does 
exactly the same job the background data collection subrou- 
tine does. It calls in subroutine to enter magnitude and 
phase data from digital voltmeters. Results are then aver- 
aged, converted to real and imaginary parts and passed to 
the main program (lines 1890 through 2090). Having both 
background and target data, main program computes system 
transfer function (lines 1040 through 1240). That concludes 
calibration procedure. The program is listed bellow. 
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10 ! “CAL IB. DRIVE©" 

20 ! 

30 ' 

40 i PROGRAM TO COMPUTE SYSTEM 
50 ' TRANSFER FUNCTION USING A 
60 ! PERFECTLY CONDUCTING 
70 ! SPHERE OVER 19-15 GHZ AT 
80 ! .1 GHZ STEPS. COMPUTATION 

90 ' BASED ON THE THEORETICAL 
190 ! VALUES OBTAINED USING THE 
110 ! PROGRAM "SPHERE. DRIVE©" . 

1 29 ' 

ISO ! THE RESULTED SYSTEM TRANS 

149 ! FER FUNCTION IS STORED AS 

150 ! 

160 ! 

170 ! "CAL I BO . DRIVE1 " C 0.875" ) 
ISO i "CALIBl . DRIVE1 " (1.125") 
190 ! "CALIB3.DRIVE1" (3. 187" > 
200 ! 

210 ! 

229 ! Cf IS THE FILE STORING 

230 ! THE SYS. TRANSF . FUNCTION 

249 ! 

250 ! Ht IS THE FILE STORING 
260 ' THEORETICAL DATA 

270 » OF THE SPHERE 

289 ! 

290 ! St DESCRIBES THE SPHERE 
3 0 0 ' 

3 l’ O C$= " C AL I B3 . DR I VE 1 " 

329 H$= " THE0R3 . DRIVE1" 

330 S $ = " 3 . 137 INCH SPHERE" 

340 A*="BKGRND. DRIVE 1 " 

350 ! 

360 X9=2*PI ! A PARAMETER 

370 ! 

380 OPTION BASE 1 
390 N9=10 ! DIMENSION OF WO 

400 F9= . 1 ! FREQ. STEP IN GHZ 

410 Ml =51 ! Ml=(U9-L9>/F9+2 

420 L9= 1 0 ! LOWER FREQ. IN GHZ 
430 U9= 1 5 ! UPPER FREQ. IN GHZ 
448 DIM A < 5 1 , 2 ) ! BACKGROUND 

OATH 

450 DIM B(51,2.) ! TARGET 

466 DIM G3(51,2) ! THEORY 

470 ! CREATE C* . Ml , 16 
480 ' CREATE A$,M1,1S 
490 ! STORE CALIBRATION DATA 
IN A FILE OF 51 RECORDS 
500 ! EACH RECORD CONTAINS ONE 
MAGNITUDE AND ONE PHASE 
510 ! DATA AT A FREQUENCY 
520 ! 

530 ASSIGN# 1 TO H* 

540 KO= < L9-2-F9 ) -T50 
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550 FOR 1=1 TO Ml 
560 K0=K0 + 50:*FS 

57© READ# 1 , K0 ; G3 Cl , 1 > , G3 C 1 , 2 > 

580 NEXT I 

59W ASSIGN# 1 TO * 

hPifl GO SUB 2160 

610 DISP "DO YOU WANT TO USE THE 
MOST RECENT BACKGROUND DATA 
? Y N " 

620 INPUT P* 

630 IF P $ = " N " THEN 680 
640 ASSIGN# 4 TO A? 

650 READ# 4 ; AC,) 

660 ASSIGN# 4 TO * 

670 GOTO 860 
680 CLEAR 

690 REMOTE 7 ! REMOTE ALL 

DEVICES 

70O CLEAR 7 ! CLEAR ALL DEVICES 
710 OUTPUT 719 "P" , L9, "21K0L3M0 

N60 1 " 

720 CLEAR 

730 DISP "REMOVE TARGET FROM CHA 
MBER , PUSH ’CONT' WHEN READY" 
740 LOCAL 7 
750 BEEP i? BEEF 
760 PAUSE 
770 REMOTE 7 
780 CLEAR 

790 DISP "TAKING BACKGROUND DATA 

li 

300 PRINT 

810 ! PRINT "BACKGROUND DATA" 

820 PRINT 

830 OUTPUT 719 s " P " , L9 , " Z 1 K0L3MO 
N60 1 " 

840 WAIT 200 
850 GOSUB 1306 
860 CLEAR 
870 LOCAL 7 

880 DISP "PUT TARGET INTO CHAMBE 
R , PUSH ’CONT' WHEN READY" 

89© DISP "TARGET IS " , S$ 

908 BEEP 0 BEEP 
910 PAUSE 
920 REMOTE 7 
930 ASSIGN# 3 TO A* 

940 PRINT# 3 ; AC,) 

950 ASSIGN# 3 TO * 

968 CLEAR 

970 DISP “COMPUTING TARGET DATA" 
980 PRINT 

990 ! PRINT "TARGET DATA" 

1000 PRINT 

1010 OUTPUT 719 ; "P" , L9, "21 K0L3M 
0N6O1 " 

1020 WAIT 200 
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1030 
i 0 *4 0 
1050 
1 060 
1 070 
1080 
1 090 
1 1 00 
1110 
1 120 
1 130 
1 140 
1 150 
1 180 
1 170 



1 180 
1 190 
1200 
1210 

1 220 
1 236 
1 240 
1250 
1260 
1270 
1 280 
1290 
1 300 

1310 

1320 

1 330 

1340 

1 350 
1 7^80 
1370 

1380 

1390 
1400 
1410 
1 420 
1 430 
1 440 
1450 
1 460 



1 470 
1480 



1 490 



GOSUB 1890 
PRINT " u 

PRINT "TRANS . FUNCT I ON " , S$ 
PRINT " " 

ASSIGN# 2 TO C$ 

FOR M=1 TO Ml 
N 1 = B ( M , 1 > - A < M , 1 > 

N 2 = B < M , 2 > - A ■' M .■ 2> 

X6=G3 < M , 1 > ✓ <: N 1 a 2+N2 a 2> 

X7 = G3 < M , 2 > -ATH2 < H2 , N 1 > 
X7=X7-X9*IHT<X7,'X9> 

IF X7> P I THEN X7=X7-X9 
PRINT# 2,M j X6 , X7 
! PRINT USING 970 ; M,X6,X7 
IMAGE DD, IX, " X6= " , SD . DDDE , 1 
X, ,, X7=", SD .DDDE 
NEXT M 

ASSIGN# 2 TO t 
CLEAR 

DISP "CALIBRATION COMPLETED 
•DATA STORED IN" , C$ 
BEEP 0 BEEP 0 BEEP 
LOCAL 7 
END 



! BACKGROUND data 
COLLECTION SUBROUTINE 

! 

! OUTPUT 10 TO 15 GHZ AT 
0 1 GHZ STEPS 

J=10*<L9-F9> ! FREQ. STARTS 

AT 10 GHZ =10040.1 GHZ 
FOR K= 1 TO Ml ! NUMBER OF 
FREQUENCY STEPS 

J=.J+ 1 0TF9 

IMAGE 1 A • 3Z , 1 4R 

OUTPUT 719 USING 1360 ; "P" 

, J, " 00Z 1 K0L3M0N6O 1 " 

! 50 MSEC NR IT FOR 

FREQUENCY TO STABILIZE 
WAIT 50 

! TAKE DATA IN FROM 722 
GOSUB 1550 
R 1 =A 1 TCQS CPI) 

I1=A1*SIN<P1> 

A<K, 1>=R1 
A < K , 2 > = 1 1 

! PRINT USING 2040 A<K,1> 

, A < K , 2 > 

NEXT K 

OUTPUT 719 ; "P" , L9, " Z1K0L3M 
0 N 6 0 1 " 

RETURN 
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1 F 0 fi ! 

1510 ! 

1520 ! 

1530 ! 

1540 ! 

1550 ! SUBROUTINE TO ENTER 

MAGNITUDE AND PHASE DATA 
FROM DIGITAL VOLTMETER 

1560 ! 

1570 ! PREPARE DIG. VOLTMETERS 
1580 ! TO SEND MAGNITUDE 
1598 ! AND PHASE DATA 
1 A 00 ! 

1610 ! N0 READINGS TAKEN AND 
AVERAGED FOR ONE ENTRY 
1620 V1=0 
1630 F 1 =0 

1640 FOR L= 1 TO N0 
1650 OUTPUT 720 " P0F 1 R 1 T 1 Z 1 FLOM 

0 “ 

1660 WAIT 10 
1670 ENTER 720 ; V 
1680 WAIT 10 

1690 OUTPUT 722 i " F 1 R7T 1 M3A0H 1 M 

1700 WAIT 10 

1710 ENTER 722 ; F 

1720 V1=V1+V 

1730 F 1 =F 1 +F 

1740 WAIT IQ 

1750 NEXT L 

1760 V=V1.'N0 

1770 F=F1/N0 

1780 A1=10 A V 

1790 P1=100TF 

1800 P1=DTR<P1> 

1810 ! PRINT USING 1310 ; K,A1,P 
1 

1820 IMAGE DD, 2X, "A= , ‘ ,MD .DDDE,2X 
, ° P= " / SD . DDDE 
1830 RETURN 
1840 ! 

1850 ! 

I860 ! 

1870 ! 

1880 ! 

1890 ! TARGET DATA COLLECTION 
SUBROUTINE 

1 900 ! 

1910 ! OUTPUT 10 TO 15 GHZ 
AT 0.1 GHZ STEPS 
1920 J=10*<L9— F9> ! FREQ . STARTS 
AT 10 GHZ =108X0.1 GHZ 
1930 FOR K= 1 TO Ml ! NUMBER OF 

FREQUENCY STEPS 

1940 J=J+10*F9 
1950 IMAGE 1A,3Z,14A 
I960 OUTPUT 719 USING 1950 ; " P" 

J. “00Z1K0L3M0N6O1" 
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1970 ' 50 MSEC WRIT FOR 

1980 ! FREQUENCY TO STRBILIZE 
1998 ' TftKE DATA IN FROM 722 
2000 GO SUB 1 550 
2010 R1=A1*C0S<P1> 

2020 I l=fil *SIN<P1 > 

2030 B<K,1>=R1 
2040 B<Ki2)=Il 

2050 ! PRINT USING 2048 ; B<K,1> 
, B < K , 2 > 

2068 IMAGE 4X, "R=" , SD . DDDE, 2X, " I 
=“ , SD . DDDE 
2078 NEXT K 

2088 OUTPUT 719 "P" , L9, "Z1K8L3M 

0 N 6 0 1 " 

2098 RETURN 

2100 ! 

2110 ! 

2120 ! 

2138 ! 

2148 ! HERDER SUBROUTINE 
2150 ! 



2150 


PRINT 


If ll 


2170 


PRINT 


II ll 


2 1 80 


CLEAR 




2190 


DISP " 
, ST 

DISP " 


CALIBRATION STANDARD" 


2280 


ENTER TGT RNG IN M" 


221G 


INPUT 


Rf 


2220 


PRINT 
" .8$ 


"CALIBRATION STANDARD 


2230 


PRINT 


"TARGET RANGE " , RT 


2240 


PRINT 


" ********** ******** " 


2259 


PRINT 


"******************“ 


22 50 


PRINT 


II II 


2270 


CLEAR 




2280 


RETURN 
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C. TARGET PROGRAM 



This program computes back scattering cross section and 
phase of the target signal return, using calibration data 
obtained in the program "CALIB . DRIVEO" and stored in the 
file C$ . Results are stored in the file G$. Comments and 
necessary program variable declarations are given at the 
begining of the main program (lines 10 through 400). 
Calibration results are then read into computer memory, as 
well as the background data (lines 410 through 550). If the 
user wants to take new background data, program calls in 
background data collection subroutine and the subroutine to 
enter magnitude and phase data from digital voltmeters. 
Target is then placed into the anechoic chamber and the 
target data collection subroutine is called along with 
subroutine to enter magnitude and phase data from digital 
voltmeters. Procedures for both background and target data 
collection are described in the calibration program section 
in detail. Having calibration, background and target data 
collected, the main program computes target cross section 
and phase shift (lines 810 through 930). Cross section and 
phase shift data could be printed, plotted or both, stored 
in a newly created or an existing file. New results could be 
averaged with the measurement results obtained and stored 
earlier. This procedure allows us to reduce system noise, 
and can be repeated as many times as we want to (lines 940 
through 2100). Header subroutine (lines 2990 through 3190) 
is written to improve data readability and to provide their 
documentation . 

MAGNITUDE PLOTTING SUBROUTINE and PHASE PLOTTING 
SUBROUTINE serve to plot target cross section and phase 
shift respectively, from the results obtained in this 
program, using HP 7225B plotter (lines 3230 through 4960). 
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On the other hand, if we want to plot earlier obtained 
and stored results, "PLOT . DRIVEO" program has to be used. 
This program is developed using above mentioned plotting 
subroutines. Target and plotting program are listed bellow. 
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10 ! "TARGET . DRIVES 11 
20 ! 

7,0 i 

4 fi i PROGRAM TO COMPUTE TARGET 

59 ! BACKSCATTERING CROSS 

60 ! SECTION AND PHASE USING 
70 ! Ct DATA AND TO STORE 

SO ! RESULTS IN THE FILE G$ . 

90 ! 

1 00 ! 

110 ! FREQUENCIES = 10-15 GHZ AT 
0.1 GHZ STEPS 

1 20 ! 

130 ! 

140 ! FILE CS STORES THE SYSTEM 
TRANSFER FUNCTION 

150 ! 

168 ! FILE G $ STORES TARGET DATA 
OBTAINED FROM THIS PROGRAM 

170 ! 

180 ! FILE H$ STORES THEORETICAL 
VALUES FOR PLOTTING OVERLAY 

190 ! 

208 C$="CALIB3. DRIVE 1 " 

210 G$= " SCR . DRIVE 1 " 

220 H$="THEGR3 . DRIVE 1 " 

230 A$= H BKGRND . DRIVE 1" 

240 ! CREATE G$,52,24 
250 ! CREATE A*, 51. -IS 
260 ! 

270 ! 

280 OPTION BASE 1 

290 N0=16 ! DIMENSION OF WO 

380 DIM R<51j2> 

310 DIM B < 5 1 .. 2 > ! TARGET 

328 DIM G4-C51..2) ! CALIBRATION 

338 DIM N ( 5 1 .. 3 ) ! RESULT 

340 DIM M9C51.3) 

350 N 1 = 5 1 ! HUMBER OF FREQ. 

CHECKED 

360 F9= . 1 ! FREQ. STEPS 

370 U9= 1 5 ! UPPER FREQ. IN GHZ 
380 L9= 1 0 ! LOWER FREQ. IN GHZ 
390 DIM T (800.. 2 > ! STORES 

THEORETICAL DATA 
400 X9=2*PI 
410 ASSIGN# 1 TO C$ 

420 READ# 1 .: G4 0> 

430 ASSIGN# 1 TO * 

440 ! MAT PRINT USING 330 ; G4 

450 IMAGE 2X,3D.4D 

460 REMOTE 7 ! REMOTE ALL 

DEVICES 

470 CLEAR 7 ! CLEAR ALL DEVICES 
480 OUTPUT 719 ; " P 1 Z 1 K0L3M0N6O 1 " 

! INITIAL SETUP OF 719 
490 CLEAR 
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500 DISP "DO YOU WANT TO USE THE 
MOST RECENT BACKGROUND 
DATA? V H " 

510 INPUT P$ 

520 IF P f = " M " THEN 576 
530 ASSIGN# 4 TO A$ 

540 READ# 4 A < , ) 

550 ASSIGN# 4 TO * 

560 GOTO 70O 

570 DISP " REMOVE TGT FROM CHAMBE 
R PUSH ’CONT’ WHEN READY" 

580 LOCAL 7 
590 BEEP 0 BEEP 
680 PAUSE 

610 DISP "TAKING BKGRND DATA" 

620 REMOTE 7 

630 OUTPUT 719 " P " , L9 , " Z 1 K0L3M0 

N60 1 " ! INITIAL SETUP OF 719 

640 WAIT 100 
65© GOSUB 2160 
660 ASSIGN# 5 TO A* 

670 PRINT# 5 A<,> 

680 ASSIGN# 5 TO t 
690 CLEAR 

700 DISP "PUT TARGET INTO CHAMBE 
R .• PUSH ‘CONT’ WHEN READY" 

710 LOCAL 7 
720 BEEP g BEEF- 
730 PAUSE 
740 REMOTE 7 

750 OUTPUT 719 ; "P" , L9, "Z1K0L3MO 
N60 1 " ! INITIAL SETUP OF 719 

760 WAIT 500 
770 GOSUB 3030 
780 CLEAR 

790 DISP "COMPUTING TARGET DATA" 

800 GOSUB 2730 

810 F0=L9-F9 

820 FOR M= 1 TO N1 

830 F0=F0+F9 

840 N(M,1>=F0 

850 X7=B<M, 1 >-A<M, 1 > 

860 X8=B<M,2>-A<M,2> 

870 X6=<X7 A 2+X8 A 2>*G4<M, 1 > 

380 N<M,2>=X6 

890 X8=ATN2<X8, X7>+G4<M, 2> 

900 X8=X3-X94: 1 NT < X8/X9 ) 

910 IF X8> P I THEN X8 = X8-X9 
920 NCM, 3 > = X 8 
930 NEXT M 

940 DISP "PRINT DATA? Y/N" 

950 BEEP 6 BEEF- 
960 INPUT PS 
970 IF P*="N" THEN 1O10 
980 PRINT " FREQ CRSEC 

PHASE" 

990 MAT PRINT USING 10O9 ; N 
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1000 

1 0 1 0 

1 028 

1 030 

1 040 

1050 

1 060 

1 070 

1080 

1090 

1 100 

1110 

1 120 

1 130 

1140 

1 150 

1 160 

1170 

1 180 

1 190 

1200 

1210 

1220 

1230 

1 240 

1 250 

1260 

1270 

1280 

1290 

1300 

1310 

1 320 

1330 

1348 

1358 

1360 

1378 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1 450 



IMAGE 2X.-3D 40 

CLEAR 

LOCAL 7 

DISP "PLOT MAGNITUDE FOR 

THIS MERSURMENT? Y/N “ 
INPUT P* 

IF P$= " H " THEN 1090 
DISP "SELECT PEN. PUSH 
' COHT ' WHEN READY" 

PAUSE 

GOSUB 3236 
CLEAR 

DISP "PLOT PHASE FOR THIS 
MEASURMENT ? Y/N" 

BEEP 6 BEEP 
INPUT Pf 

IF PT= " H " THEN 1170 
DISP "SELECT PEN. PUSH 
* CONT' WHEN READY" 

PAUSE 

GOSUB 4340 
CLEAR 

DISP "DO YOU WANT TO MAKE 
AVER. WITH PREVIOUS DATA?" 
DISP "?Y/N" 

BEEP 0 BEEP 
INPUT P* 

IF P$=“Y" THEN 1426 
DISP "DO YOU WANT TO STORE 
DATA ? Y/N " 

INPUT P* 

IF P$="N" THEN 1886 
M0= 1 

DISP "DO YOU WANT TO STORE 
DATA IN FILE" 

DISP G$ 

DISP "? Y/N" 

INPUT P* 

IF P$="Y" THEN 1388 
DISP "ENTER NAME OF THE DAT 
A FILE TO BE USED FOR STORA 
GE" 

INPUT G* 

DISP “IS THIS AN OLD FILE 
TO BE UPDATED ? Y/N " 
INPUT P* 

IF P$="Y° THEN 1386 
CREATE GT, 52,24 
ASSIGN# 2 TO G* 

PRINT# 2 i M0,N-w) 

ASSIGN# 2 TO * 

GOTO 1880 

DISP “ARE THE DATA STORED 
IN FILE" 

DISP G* 

DISP "? Y/N" 

INPUT P$ 
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1460 IF P$="Y" THEN 1490 
1470 CUSP "ENTER NAME OF DATA 
FILE TO EE USED FOR 
THE AVERAGE" 

1480 INPUT GT 

1496 ASSIGN# 6 TO G$ 

1506 READ# 6 MO, M9< , > 

1510 ASSIGN# 6 TO * 

1520 FOR K= 1 TO N1 
1 530 M9 < K , 2 > = M9 < K , 2> TM0+N < K , 2 > 
1540 M9 ■: K 2 > = M9 < K , 2 ) ✓ C M0+ 1 > 

1 550 M9 < K , 3 > = M9 C K , 3 > TMO+N < K .. 3 > 

1 560 M9 C K , 3 > = M9 < K , 3 > / < MO+ 1 > 

1570 NCK, 2>=M9<K, 2> 

1580 N < K, 3 ) =M9 ( K , 3 > 

1590 NEXT K 

1600 N0=M0+1 

1610 ASSIGN# 7 TO G$ 

1620 PRINT# 7 ; M 0 , N < , ) 

1638 ASSIGN# 7 TO * 

1640 PRINT "DATA IS AVERAGE OF", 
M0, "MEASURMENTS" 

1650 DISP "PRINT DATA? Y-'N" 

1668 BEEP S BEEP 
1670 INPUT Ft 
1680 IF P $ = " N " THEN 1720 
1690 PRINT " FREQ CRSEC 

PHASE" 

1700 NAT PRINT USING 1000 ; N 
1710 INAGE 2 X , 3 D . 4 0 
1728 CUSP "PLOT MAGNITUDE? Y/'N" 
1730 BEEP 6 BEEP 
1740 INPUT Ft 
1750 IF P $ = " N " THEN 1790 
1760 DISP "SELECT PEN FOR MAGNIT 
UDE PLOT. PUSH * CONT ■ WHEN 
READY . " 

1770 PAUSE 
1780 GOSUB 3230 
1790 CLEAR 

1800 CUSP "PLOT PHASE? Y/N" 

1810 BEEP 6 BEEP 
1820 INPUT Ft 
1830 IF P*="N" THEN 1870 
1840 DISP "SELECT PEN AND CHANGE 
PAPER FOR PHASE PLOT. 
PUSH 'CONT 1 WHEN READY." 
1850 PAUSE 
I860 GOSUB 4340 
1870 CLEAR 

1880 DISP "DO YOU WANT TO" 

1890 DISP "OBTAIN DATA" 

1900 DISP "FOR A NEW TARGET?" 
1910 DISP " " 

1920 CUSP "ENTER Y.-'N" 

1930 INPUT Ft 

1940 IF P $ = " N " THEN 2070 
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OISP "DO YOU WANT TO USE 
THE SAME FILE" 

DISP Gf 

DISP "TO STORE NEW ORTA? 
Y-'N" 

INPUT P* 

IF P$="Y" THEN 2066 
DISP "ENTER NEW FILE NAME 
TO STORE TARGET DATA" 
INPUT Gt 

DISP "IS THIS AN OLD FILE 
TO BE UPDATED? YVN" 
INPUT P* 

IF p$="Y" THEN 2066 
CREATE G$ / 53 > 24 
GOTO 460 
CLEAR 

DISP "END OF PROGRAM " 

BEEP A BEEP e BEEP 
END 



! BACKGROUND DATA 
! COLLECTION SUBROUTINE 
! 

! OUTPUT ID TO 15 GHZ 
J=10*<L9-F9> ! FREQ. STARTS 

AT 10 GHZ TO BE INCREM. 
AT 0.1 GHZ STEPS 



FOR K= 1 TO N1 ! NUMBER OF 
FREQUENCY STEPS 



J= J+ 1 0$F9 

IMAGE lft,3Z,14A 

OUTPUT 719 USING 2228 ; "P" 

, J, "00Z1 K0L3M0N6O1 " 

! 50 MSEC WAIT FOR 

FREQUENCY TO STABILIZE 
WAIT 50 

! TAKE DATA IN FROM 722 
GOSUB 2420 
R1=A1*C0S<P1 ) 

Il=fll*SIN<Pl> 

A<K, 1 ) =R 1 
ft < K , 2 > = 1 1 

! PRINT "I1=",A(K,2> 

! PRINT "Rl=",fi<K, 1) 

NEXT K 

OUTPUT 719 "P" , L9, "Z1K0L3M 

0N6O 1 u ! IHIT. SETUP OF 719 
RETURN 
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2 4 0 y 
2 4 1 0 
2420 



2430 

2440 



2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2680 

2610 

2628 

2630 

2640 



! SUBROUT I HE TO ENTER MAGNI 
TUOE AND PHASE DATA FROM 
DIGITAL VOLTMETERS 

i 

! PREPARE DIGIT. VOLTMETERS 
TO SEND DATA 
! H0 READINGS TAKEN AND 
AVERAGED FOR ONE ENTRY 
V1=0 
wi=0 

FOR L= 1 TO H0 

OUTPUT 720 ; "F1R1T1Z1 FL0MO" 
WAIT 10 

ENTER 720 ; V6 
WAIT 18 

OUTPUT 722 ; "F1R7T1 M3A0H1 " 
WAIT 10 

ENTER 722 ; W0 

V 1 =V 1 +VO 

Wl=Wl+WO 

WAIT 10 

NEXT L 

V0=V1/N8 

W0=W1/H0 

A 1 = 1 0 A V0 

P 1 = 1 00$W0 

! PRINT " A 1 = *' , A 1 



2650 PI =DTR <P 1 > 

2660 ! PRINT U P1 = %P1 
2670 RETURN 
2680 ! 

2690 ! 

2708 ! 

2710 ! 

2720 ! TARGET DATA 

2730 ! COLLECTION SUBROUTINE 
2740 ! 

2750 ! 

2760 ! OUTPUT 10 TO 15 GHZ AT 
0.1 GHZ STEPS 

2770 -J=10*<L9-F9> ! INITIAL FREQ 

UENCY AT 10 GHZ=100*0.1 GHZ 
2780 FOR K= 1 TO N1 ! FREQ. STEPS 
2790 J=J+10*F9 ! .1 GHZ IHCREM. 

2890 IMAGE 1A,3Z,14A 
2810 OUTPUT 719 USING 28O0 .; "P" 

, J, u 0OZ 1 K0L3M0N60 1 " 

2820 ! 50 MSEC WAIT FOR 

FREQUENCY TO STABILIZE 
2830 WAIT 50 

2840 ! TAKE DATA IN FROM 722 
2850 GOSUB 2428 
2860 R 1 =A 1 YCOS < P 1 ) 

2870 I 1=A1*SIN<P1> 

2880 B<K,1)=R1 
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2*96 B < K , 2 > = I 1 

2*00 ! PRINT " R 1 = " .• B v K .• 1 ) 

2*10 ! PRINT " 1 1 = ", BCK, 2) 

2*20 NEXT K 

oc.70 riUTPUT 71* : “P" , L9, "21K0L3M 
0N601 '• ! I HIT. SETUP OF 71* 

2*40 RETURN 
2*50 ! 

2*60 ! 

2970 ! 

2*80 ! 

2**0 ! HERDER SUBROUTINE 
3000 ! 

3010 ! 

3020 D * = " M 0 H T H / D A T E / Y E A R " 

3030 PRINT " " 

3040 PRINT " " 

3050 BLEAR 

3060 DISP "ENTER TODAY'S DATE - 
MONTH, DATE, YEAR" 

3070 INPUT D* 

3080 DISP "ENTER TGT DESCRIPTION 

ll 

30*0 INPUT T* 

3100 DISP "ENTER TGT RANGE" 

3110 INPUT R$ 

3120 PRINT DT 

3130 PRINT "TARGET IS " , T$ 

3140 PRINT "TGT RANGE", R$ 

3150 PRINT "******************" 
3160 PRINT "******************" 
3170 PRINT " " 

3180 CLEAR 
3190 RETURN 
3200 ! 

3210 ! 

3220 ! 

3230 ! MAGNITUDE PLOTTING 
3240 ! SUBROUTINE 
3250 ! 

3260 ! 

3270 PLOTTER IS 705 
3280 LOCATE 32,122,20,85 
3290 FRAME 

3300 ! SEARCH FOR MAX Sc MIN. 
3310 ! S0=N <1,2 > 

3320 ! S1=S0 
3330 ! FOR M=2 TO Ml 
3340 ! IF S0> N ( M , 2 > THEN S0=N<M, 
2 > 

3350 ! IF SI <N<M, 2> THEN S1=N<M, 
2 ) 

3360 ! NEXT M 

3370 DISP "ENTER LOWER VALUE FOR 
MAGNITUDE PLOTTING" 
3380 BEEP 6 BEEP 
3390 INPUT S0 
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3400 


DISP "ENTER UPPER VALUE 


FOR 




MAGNITUDE PLOTTING 


n 


3410 


INPUT SI 




3420 


L 1 = I NT L9> 




3430 


IJ1=CEILCU9> 




3440 


S3=LGT < S 1 > 




3450 


S4= I NT < S3 > - 1 




3460 


.'35=S3 — 84 




3470 


85= I NT < 1 0^85 >+l 




3430 


84=10^84 




3490 


L0=INT <S0'S4) 




3500 


IF 85— L0<=1 4 THEN 3580 




3510 


IF S5-L0>=50 THEN 3550 




3520 


85= . 5*85 




3530 


84=841:2 




3540 


GOTO 3498 




3550 


85= . 2*85 




3560 


84=5*84 




3570 


GOTO 3490 




3530 


L0=S4*L0 




3590 


IJ0 = S5*S4 




3600 


D0=U0-L0 




3610 


SCALE L1,U1,L0,IJ0 




3620 


FXD 0 , 4 




3630 


LAKES -1,S4,L1,L0 




3640 


MOVE LI.. 0 




3656 


FOR K= 1 TO HI 




3660 


M9=N < K , 1> 




3670 


R 9 = N < K , 2> 




3630 


MOVE M3 , R9 




3690 


DRAW M9,R9 




3700 


! GOSUB 4060 




3710 


NEXT K 




3720 


M5=<U1+L1 >-'2 




3730 


MOVE M5,L0-.09*D0 




3740 


LORG 5 0 CSIZE 3,1,0 




3750 


! LABEL "FREQUENCY <GHz 


> " 


3760 


MOVE M5,L0-.16*D0 




3770 


! LABEL "Fig. . . TARGET2 




- Cross Section" 




37S0 


MOVE L1-.7, .5*CL0+U0> 




3790 


LDIR PI, -2 




3800 


! LABEL "CROSS SECTION 
METERS >" 


<SQ. 


3810 


MOVE M5, IJ0+ . 895*00 




3820 


LDIR 0 




3ft30 


CSIZE 3,1,0 




3840 


! LABEL "TARGET: SPHERE 
.125 inch." 


D=1 


3850 


MOVE M5,IJ0+.06*D0 




3860 


! LABEL "POLARIZATION: 
ZONTAL" 


HOP. I 


3870 


MOVE M5,U0+ .025*00 




3830 


! LABEL "MATERIAL: ALUMINUM 

II 


3890 


PEN UP 
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3906 DISP "OVERLAY THEORETICAL 
CURVE? Y/N" 

3910 ! 

3920 INPUT Pf 

3930 IF Pi=" N " THEN 4190 

3940 DISP "IS THE THEORETICAL 

DATA STORED IN THE FILE" , H$ 
3950 DISP "? Y/N" 

3960 INPUT P$ 

3970 IF P$=" Y" THEN 4000 
3980 DISP "ENTER NAME OF THE 

DATA FILE TO BE PLOTTED." 
3990 INPUT HT 
4000 BEEP 0 BEEP 

4010 DISP "CHANGE PEN IF DESIRED 
. PUSH 'CONT' WHEN READY." 
4020 PAUSE 

4030 ASSIGN# 3 TO H$ 

4040 Jl=<L9-2>*56 

4050 J2=CU9-2>*50 

4060 FOR J=J1 TO J2 

4070 READ# 3,J ; T< J , 1 > , T < J , 2 > 

40-80 NEXT J 

4090 ASSIGN# 3 TO X 

4100 F0=L9 

4110 P.9 = T < J 1 , 1 > 

4120 HOVE F0,R9 
4130 FOR I=J1 TO J2 
4140 F0=F0+ . 02 
4150 R9=T<I, 1> 

4160 DRAW FO , F‘9 
4170 NEXT I 
4180 PENUP 
4190 RETURN 
4200 ! 

4210 ' 

4220 ! 

4230 ! PLOT CROSS 
4240 MOVE M9,R9 
4250 CSIZE 2, .5,0 
4260 LABEL "+" 

4270 ! I MOVE .00025- 00825 
4280 ! I DRAW -.0005,0 
4290 RETURN 
4300 ! 

4310 ! 

4320 ! 

4330 ! 

4340 ! PHASE PLOTTING SUBROUTINE 
4350 ! 

4360 PLOTTER IS 705 

4370 LOCATE 32,122,20,85 

4380 FRAME 

4390 S4= . 5 

4400 110=1 

4410 L0= - 1 

4420 D0=U0-L0 
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4430 SCALE L 1 , U 1 , L0 , U0 
4440 F X D 0,3 
4450 LAXES -1,S4,L1,L0 
4460 MOVE LI , 0 
4470 FOR K= 1 TO HI 
44S0 M9=N<K,1> 

4490 R9 = N<K,3>,'PI 
4500 MOVE M9,P.9 
4510 DRAW M9,R9 
4520 ! GOSUB 4680 
4530 NEXT K 

4540 MOVE M5, L0- . 09*D0 
4550 LORG 5 & CSIZE 3,1,0 
4560 ! LABEL "FREQUENCY <GHz>" 
4570 MOVE M5, L0-. 16TD0 
4580 ! LABEL “Fi<?. . . TARGET2 

- Phase" 

4590 MOVE L1-.7, .5*<L0+UCO 

4600 LDIR PI-' 2 

4610 LABEL "PHASE <PI>" 

4620 MOVE M5,U0+.095*D0 
4630 LDIR 0 
4640 CSIZE 3,1,0 

4650 ! LABEL "TARGET = SPHERE D=1 
.125 inch." 

4660 MOVE M5..IJ0+ 06*D0 
4670 ' LABEL "POLARIZATION: HORI 
ZOHTAL " 

4680 MOVE M5, U0+ . 025TD0 

4690 ! LABEL "MATERIAL: ALUMINUM 

II 

4700 PENUP 

4710 DISP "OVERLAY THEORETICAL 
CURVE? Y.'N" 

4720 INPUT PT 

4730 IF P*="N" THEN 4860 

4740 BEEP 0 BEEP 

4750 DISP "CHANGE PEN IF DESIRED 
. PUSH ’ CONT ‘ WHEN READY." 
4760 PAUSE 
4770 F0=L9 
4780 R9=TC J1 , 2>.-'PI 
4790 MOVE F0,R9 
4800 FOR I=J1 TO J2 
4810 F0=FA+ 02 
4828 R9=T <1,2 >/PI 
4830 DRAW F0,R9 
4840 NEXT I 
4850 PENUP 
4860 RETURN 
4870 ! 

4880 ! 

4890 ! 

4900 ! PLOT DOT 
4910 MOVE M9,R9 
4920 CSIZE 2, .5,0 
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493G LRBEL 

4940 ! I MOVE . 00025 , . 00025 
4950 ! I DRAW -.0005,0 
4960 RETURN 
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! "PLOT . DRIVED" 

i 

' PROGRAM TO PLOT 
! BACKSCATTERING CROSS 
! SECTION FIND PHASE VS. 

! FREQUENCY USING 
! STORED DATA. 

| 

G$=“SCR4 DRIVE1 " 

H$="THEOR0. DRIVE 1" 

OPTION BASE 1 
DIM N < 5 1 , 3 :> ! DATA TO BE 

PLOTTED 

N 1=51 ! NUMBER OF 

FREQUENCES CHECKED. 

F3= . 1 ! FREQ. STEPS 

IJ9=1 5 ! UPPER FREQ. IN GHZ 
L3= 1 0 ! LOWER FREQ. IN GHZ 
DIM T< 809,2) ! STORES 

THEORETICAL DATA 

CLEAR 

DISP "DO YOU WANT TO PLOT" 
DISP "FROM THE F I LE : " 

DISP G$ 

DISP "ENTER Y/N?" 

INPUT PT 

IF P $ = " Y " THEN 230 

DISP "ENTER NAME OF FILE “ 

DISP "TO BE PLOTTED" 

BEEP 0 BEEP 
INPUT Gf 
ASSIGN# 1 TO G$ 

READ# 1 ; M0 , N C , ) 

ASSIGN# 1 TO * 

X9=2*PI 

DISP "PRINT DATA? Y/N“ 

BEEP 0 BEEF- 
INPUT P$ 

IF P f = " H " THEN 460 
PRINT " FREQ CRSEC 

PHASE" 

MAT PRINT USING 330 .; N 

IMAGE 2X, 3D . 4D 

CLEAR 

DISP "PLOT MAGNITUDES Y/N" 
BEEP 0 BEEF- 
INPUT P* 

IF P $ = " N " THEN 430 
DISP "SELECT PEN FOR MAGNITU 
DE PLOT. PUSH ' CONT' WHEN RE 
ADY . “ 

PAUSE 
GOSUB 740 
CLEAR 

DISP "PLOT PHASE? Y/N" 

BEEF' 0 BEEF- 
INPUT P$ 
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IF P$=" N " THEN 560 
DISP "SELECT PEN RND CHANGE 
PAPER FOR PHASE PLOT. PUSH 
1 CGNT ' WHEN READY." 

PAUSE 

GOSUB 13O0 
CLEAR 

DISF "DO YOU WANT TO" 

DISP "PLOT DATA" 

DISP "FORM ANOTHER FILE?" 
DISP " " 

DISP "ENTER Y/N" 

INPUT FT 

IF P$=“Y" THEN 110 
CLEAR 

DISF "END OF PROGRAM" 

BEEP 6 BEEP & BEEF- 
END 



! MAGNITUDE PLOTTING 
! SUBROUTINE 

i 

PLOTTER IS 705 
LOCATE 32,122,20,35 
FRAME 

! SEARCH FOR MAX. & MIN. 

! S 0 = N < 1 , 2 > 

! 31=30 

! FOR M=2 TO N1 
! IF S0>N<:M, 2) THEN S0=N<M, 
> 

! IF S 1 < N ( M , 2 > THEN 81=N<M, 

> 

! NEXT M 

DISP "ENTER LOWER VALUE FOR 
MAGNITUDE PLOTTING" 

BEEP S BEEP 
INPUT SO 

DISP "ENTER UPPER VALUE FOR 
MAGNITUDE PLOTTING" 
INPUT SI 
L 1=1 NT (LSO 
U1=CEILCU9> 

S3=LGT (SI > 

S4=INT <S3>-1 

S5=S3-S4 

S5= I NT < 10-35) + 1 

S4= 1 0-S4 

L0=INT <S0/'S4> 

0 IF S5-L0 <= 1 4 THEN 1030 
3 IF S5-L0> =50 THEN 1050 
0 S5= . 5*S5 
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163 0 
1 040 
1050 




1060 ^4“ y 4 

1070 GOTO 990 

1080 L0=S4*L0 

1090 !J0=S5TS4 

1 100 D0=U0-L8 

1110 SCALE LI , U1,L0,U0 

1120 FXD 0,4 

1130 LAKES -1,S4,L1,L0 

1140 MOVE L 1 , 0 

1150 FOR K=1 TO HI 

1 160 M9=N < K , 1> 

1170 R9=N < K , 2> 

1130 MOVE M9,R9 

il96 DRAW M9,R9 

1200 NEXT K 

1210 M5 = aJl+Ll >.'-2 

1 220 MOVE M5 , L8- . 09:4:00 

1230 LORG 5 6 CSIZE 3,1,0 

1240 LABEL "FREQUENCY (GHz)“ 

1250 MOVE M5,L0-.16*D0 
1260 LABEL "Fig. 3. 13. TARGET3 - 
CROSS SECTION" 

1270 MOVE L1-.7, 5*<L0+U0> 

1286 LOIR PI. '2 

1290 LABEL "CROSS SECTION <SQ . ME 
TEES > " 

1 308 MOVE M5, U0+ . 03TD0 
1310 LOIR 0 
1320 CSIZE 3,1,8 

1330 LABEL "TARGET : SPHERE <0=3. 

187 in . 'ALUMINUM V* 

1340 MOVE M5,U0+. 83*06 
1350 LABEL " POLAR I ZAT I ON : HORIZO 
NTAL " 

1360 PENUP 

1370 DISP "OVERLAY THEORETICAL 

rt mi irri v , 1 ^ l i *i 



1380 ! 

1390 INPUT P* 

I486 IF P*="N" THEN 1666 
1410 DISP "IS THE THEORETICAL 

DATA STORED IN THE FILE" ,H* 
1420 DISP "? Y,'N" 

1430 INPUT PS 

1440 IF P$="Y“ THEN 1470 

1450 DISP "ENTER NAME OF THE 

DATA FILE TO BE PLOTTED." 
1460 INPUT HT 
1470 BEEP 0 BEEP 

1480 DISP "CHANGE PEN IF DESIRED 
. PUSH 'CONT' WHEN READY." 
1490 PAUSE 

15O0 ASSIGN# 3 TO H$ 

1510 J1 = <L9-2V*50 



CURVE? Y N 
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1S-F! l2=''JJ9-2>'t'58 

1538 FOR J= J 1 TO J2 

1 546 READ# 3 , J ■ T < J , 1 > .■ T < J 2> 

1550 NEXT J 

1560 ASSIGN# 3 TO % 

1576 F6=L9 
1586 R9=T<J1,1> 

1596 MOVE F 6 , R9 
1666 FOR I=.J1 + 1 TO J2 
1616 F6=F6+ . 62 
1620 R9=T ■' I 1 > 

1636 DRAW F6,R9 
1646 NEXT I 
1656 PENUP 
1666 RETURN 
1670 ! 

1 686 ! 

1 696 ! 

1706 ! PLOT CROSS 
1716 MOVE M9,R9 
1720 CSIZE 2 , 5,G 
1736 LABEL " + " 

1740 ! I MOVE . 06025, . 00625 
1756 ! I DRAW -.0605,0 
1766 RETURN 
1776 ! 

1 780 ! 

1796 ! 

1806 ! PHASE PLOTTING SUBROUTINE 
1816 ! 

1820 PLOTTER IS 765 
1836 LOCATE 32,122,26,85 
1846 FRAME 
1856 34= . 5 
I860 L 1 = I NT < L9) 

1870 U1 =CEIL<U9) 

1886 U0= 1 

1896 L 6 = - 1 

I960 O0=IJ0-L0 

1910 M5= < U 1 +L 1 > s2 

1926 SCALE L1,IJ1,L0,U0 

1936 FXD 0,2 

1946 LAXES - 1 , S 4 , L 1 , L 0 

1956 MOVE L 1 , 0 

I960 FOR K= 1 TO N1 

1976 M9=N <K , 1 > 

1986 R9=N<K,3>/PI 
1996 MOVE M9,P.9 
2600 DRAW M9,R9 
2016 NEXT K 

2020 MOVE M5, L6- . 69*06 
2630 LORG 5 <*! CSIZE 3,1,6 
2640 LABEL "FREQUENCY <GHz>" 

2650 MOVE M5, L6- . 16TD0 
2960 LABEL "Fig. 3. 14. TARGET3 - 
PHASE" 

2076 MOVE L1-.7, .5*<L0+UO> 
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2686 LG IF: PI- -2 

2026 LfiBEL " PHftSE CPI)" 

2106 MOVE M5,U6+. 68*00 
2116 LDIR 6 
2126 CSIZE 3,1,0 

2136 LfiBEL "TARGET : SPHERE <D=3. 

1 87 i n . ,'fiLUM I HUM > " 

2140 MOVE M5, U0+ . 83*D8 
2150 LfiBEL " POLfiR IZfiT I OH : HORIZO 
MTfiL " 

2166 PEN UP 

2170 OISP "OVERLAY THEORET I CfiL C 
URVE? Y/N 1 ' 

2186 INPUT P$ 

2190 IF P$ = " N " THEN 2320 
2206 BEEP 6 BEEP 

2210 DISP "CHANGE PEN IF DESIRED 
. PUSH 1 CCiNT' WHEN REfiDY . " 
2228 PAUSE 
2236 F6=L9 
2246 R9=T < J 1 , 2 > - / P I 
2250 MOVE F0,R? 

2266 FOR I=J1+1 TO J2 
2270 F 6 = F 6 + . 8 2 
2286 R9=T< 1 , 2 > . y P I 
2296 DRAW F0,R9 
2306 NEXT I 
2310 PEHUP 
2326 RETURN 
2336 ! 

2346 ! 

2350 ! 

2360 ! PLOT DOT 
2376 MOVE M9,R9 
2386 CSIZE 2, . 5 , Q 
2390 LABEL " t " 

2406 ! IMOVE . 00625, . 00625 
2416 ! IDRFIW -.0005,6 
2420 RETURN 
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V. SUMMARY AND CONCLUSIONS 



Demands on target cross section measurements and anal- 
ysis have increased considerably in recent years. There are 
several reasons for this : 

• Great attention has been paid to cross section reduc- 
tion through the applications of radio frequency 
absorbing materials. Returns from many targets of prac- 
tical interest are suppressed. 

• Introduction of nonmetallic materials into target 
construction has further complicated target scattering 
properties . 

• Sophisticated radar systems, on the other hand, 
provided additional means for target discrimination. 
For example, the polarization characteristics of a 
target can be acquired and utilized. 

In this work the calibration of the CW scattering range 
was performed. The theoretical solution had been tested 

against experimental results for several perfectly 
conducting spheres of different diameters. Back scattered 
fields were measured for three spheres and five finite 
length, thin walled, right circular cylinders. 
Corresponding back scattering cross sections and return 
phase shifts were obtained. Differences between theory and 
experiment were discussed in terms of noises introduced by 
such sources as background reflections, system noise, 
target - antenna distance , and target support problems. 
Agreement between theory and measurement data was good if 
the calibration and target data acquisition were performed 
in a short time interval. Otherwise, as the time elapses, 
differences between theoretical and measurement data became 
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more and more pronounced. This problem requires further 
study . 

Measurements on the cylinders showed that for broad side 
incidence, the number and location of the resonance frequen- 
cies are strongly dependent on the diameter of the cylinder. 
The back scattering cross section depends on both the 
cylinder diameter and its length. The phase shift, on the 
other hand, depends only on the cylinder diameter. A change 
in polarization provides additional cross section and phase 
data. Several cylinders (targets) might be distinguished by 
comparing acquired data of their phases and cross sections 
for two orthogonal polarizations. Future works will inves- 
tigate applicability and reliability of this approach towad 
the identification of targets of complex shapes. 
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